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Abstract 

The renormalization evolution of all parameters in the neutrino mass matrix depends only 
on one variable, the energy scale. This fact, coupled with rephasing considerations, leads to 
a set of renormalization invariants, correlating the evolution of physical parameters. For the 
general three flavor case, we obtain these invariants explicitly and discuss their implications. 

1 Introduction 

Recent results from the atmospheric and solar experiments have shown strong evidence for 
neutrino oscillations. |2| These observations indicate that neutrinos are not massless, and 
that two of the neutrino mixing angles are large, or even maximal. This is in contrast to the 
quark mixing angles, which are all small. Even with our limited knowledge, it seems clear that 
the neutrino mass matrix, just like their quark counterpart, has a rich structure, and it is urgent 
to have an understanding of its salient features. To account for the minuscule neutrino masses, 
the seesaw model ||] makes use of a heavy scale for the right-handed neutrinos. Thus, any 
theoretical understanding of the observed neutrino parameters necessarily involves two vastly 
different energy scales. This means that renormalization effects must be taken into account in 
any theoretical model of the neutrino mass matrix. 

The renormalization of the neutrino mass matrix has been extensively discussed in the 
literature [|, |6|. Although different models, such as the SM or MSSM, give rise to numerically 
distinct results for individual parameters, as we have shown for the two flavor problem ||, there 
are renormalization group equation (RGE) invariants, correlating the evolution of the physical 
parameters. These invariants are the consequences of the general structure of the RGE, and 
remain the same for any model which conserves FCNC. 
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In general, RGE evolution implies that each physical parameter becomes a function of the 
energy scale. Thus, if there are n independent parameters in the mass matrix, we might expect to 
have (n — 1) RGE invariants. However, the mass matrices are also subject to arbitrary rephasing 
transformations. In addition, such phases are generated by the renormalization transformation. 
As a result, the physical RGE invariants must also be rephasing invariants. For the three flavor 
mass matrix, it turns out that there are three (complex) RGE and rephasing invariants, amongst 
its eight physical parameters. These invariants will be detailed in Sec. 4. Just as for the two 
flavor problem, these invariants are universal, valid for a wide class of models. 

2 Parameterization of Mass Matrices 



Before we embark on a detailed discussion of renormalization, it is useful to introduce a general 
parameterization of the mass matrix which will facilitate the analysis. In this paper, we will 
consider the (symmetric) neutrino mass matrix to originate from a dimension five term in the 
effective Lagrangian, 

L = fu T u{<j>)(<j>) = u T My (1) 

Here, ((f)) is the vacuum expectation value of the Higgs scalar, / is the coupling constant, M® 
is the neutrino mass matrix, and v is the neutrino wave function in the flavor basis. 
We can write, in general, 

M° = UM dias U T , (2) 
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Here, the mass eigenvalues are given by exp(2r/j), cc, are the unphysical phases from the neutrino 
wave functions, (e2, £5, £7) are the physical neutrino mixing angles, £3 is a CP violating phase, the 
7's are the intrinsic CP phases of the mass eigenvalues, and the A's are the Gell-Mann matrices. 
To preserve the symmetry of the flavors, we will not use the diagonal A's in M diag . Note also 
that, since neutrino oscillations are governed by the effective Hamiltonian, H = (M v Ml)/2E, 
they are independent of the phases 7,. 
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It is convenient to factor out the determinant (the overall scale) of M® and define 

M u = (detM°)- 1/3 M°, 
delM v = 1. 

The condition detM u = 1 is obtained by imposing the following relations on M®: 

Yj Vi = a i = Yl Ti = °- 



(7) 
(8) 
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This corresponds to the fact that, with detM u = 1, M v only depends on the mass ratios (A77) 
and the relative phases (Aa and A7). Note that, with the parametrization in Eq.(||), M v can be 
analytically continued into an SU(3) matrix (rjj — ► ir/j). A fact which will be used later. Also, 
rephasing of the neutrino wave functions changes only ctj, while leaving all other parameters 
invariant. 

It is also useful to write down the symmetric matrix M u explicitly. 
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Here, we use the notation S( 2 ) = sine2,s 2 (2) = sin2e2, etc. We have also set a« = 0, without 
loss of generality. 



3 RGE 



The RGE for the effective neutrino mass matrix, M u , has been studied extensively. In the SM 
and MSSM, the RGE were obtained explicitly and can be written in the form @, |5|, § 

±M° = kM° + {q,M°}, (16) 
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where k is a constant, Q is a diagonal and traceless matrix 



trQ = 0. (17) 

t is the scale variable t = ln(////io), with (/U,/xo) = energy scale. The solution is given by 

detM°(t) = e 3Kt detM°(0), (18) 
Af„(t) = e Q1 M,(0)e Qi . (19) 

The quantities k and Q were given explicitly in terms of the leptonic Yukawa constants. 

We note that the form of Eq. ( |T9| ) is quite general. The operator e® 1 originates from the (rela- 
tive) neutrino wave function renormalization. It amounts to a change of relative scale (rescaling) 
between the different flavors. As long as the interactions responsible for the renormalization 
do not contain FCNC, their effect will be a diagonal matrix multiplication, as in Eq.(|i"S|). The 
close relation between rescaling (renormalization) and rephasing is revealed by considering pure 
imaginary values for Q, which turns Eq.(^) into a rephasing transformation. If, in addition, 
we consider r\j — > ir\j in M v , then the equation becomes a rephasing transformation in SU{3). 



Eq.(19) is a formal solution of RGE, since it only gives the t-dependence of the matrix 
elements of M v . One would really like to know the i-dependence of the physical parameters. 
To this end, we must reexpress M v in Eq.(19) in the form of Eq.(S), 



M„(t) = U{t)M di ^{t)U T {t), (20) 

and one need to relate the physical parameters at scale t to those at t = 0. 

Mathematically, the RGE solution corresponds to the relation connecting the parameters, 
known as the Baker-Campbell-Hausdorff (BCH) formula, between different rearrangements (e.g. 
Eqs.(|l9|) and (p0|)) of the non-commuting factors in an element of SL{3, C). Since we are dealing 
with only exponential functions which are free of singularities, these relations should remain valid 
under analytic continuations. In particular, the BCH formulae for SL(3, C) are the analytic 
continuation of those for SU(3). These ideas can be implemented explicitly for the case of two 
flavors, which we will study first before we take on the full analysis of the three flavor problem. 

Consider a general symmetric SU(2) matrix (with two real parameters), which can be written 

as 

=e -t/9oa e 2i7os e i/9o! 2j (21) 



or 



N2 = e i " a3 e 2ifai e i " a3 . (22) 
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The relations (BCH formula) between the two parametrizations Ni and N2 can be read off from 
the matrix elements as given in Eqs (|2l"|), (|22"|). 

cos 27 = cos 2f cos 2a), (23) 

sin 27 cos 2/? = cos2fsin2cD, (24) 

sin 27 sin 2(3 = sin2f. (25) 

When N is subject to a rephasing transformation, 

N ^ e - i&<T3 Ne- i&as , (26) 



it is obvious that parametrization Eq,(^2|) makes it trivial, resulting in uj — > a) — a. But for N\ 
in Eq. (|2l|) , it induces the change (3 — > (3' , 7 — > 7', satisfying 



cos 27' = cos 2f cos 2(a) — a), (27) 

sin 27 cos 2/3 = cos 2f sin 2(a) — a), (28) 

sin 27 sin 2/?' = sin2f. (29) 

It follows that 

tan 2/3' = f Wl C ° S 2 " . (30) 

cos 2p — tan 2a/ tan 27 

In addition, f is invariant under rephasing. Thus, we have the rephasing invariant, in terms of 
the parametrization in Eq.([H]). 

sin 27 sin 2(3 = sin 27 sin 2(3'. (31) 



Eqs. (|30| ) and ( pl|) are the solutions to the SU(2) rephasing transformation, Eq.(p6|). 

The same results can be taken over for symmetric mass matrices, where all variables (/3, 7, a), f) 
are complex, with four real parameters. Let us first consider the case of real mass matrices, 
corresponding to pure imaginary (7,d>,f) : 

(7, a), f) — > (i7, iu), it). (32) 

The resulting mass matrix can be written in two alternative forms, 

N 2 = e-^e-^e-^ 3 . (33) 

With a — » ia, the rephasing transformation on N becomes a renormalization (rescaling) trans- 
formation on N : 

N — > e aa ' A Ne aa K (34) 
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The solution to the RGE is obtained directly from the BCH formula of SU(2), Eg. fl30|) , by 
analytic continuation: 



tan 2(3' 



sin 2(3 j cosh 2a 



(35) 



cos 2(5 — tanh 2a/ tanh 27 
This is the same relation obtained earlier for RGE evolution. At the same time, we have an 
RGE invariant @ : 

sin 2(3 sinh 27 = sin 2(3' sinh 27' . (36) 

Note that this is also automatically invariant under (relative) rephasing on the mass matrix N. 
In addition, while the value a in Eq.(|34]) is model dependent, Eq.(|36|) is not. 

In general, all parameters ((3,j,uj,f) are complex, so that Ni and N 2 become complex mass 
matrices. As was shown in Ref.jT]], we can demand that (3 and (3 1 be real by adding rephasing 
factors to N\. The resulting generalization of Eq.(|30|) coincides with Eq.(15) of Ref.||. In this 
connection, we emphasize that the renormalization transformation, Eq.(|34|), generates rephasing 
factors when N is complex, as was shown explicitly in Ref.||. 

Schematically, solutions to the rescaling transformations on mass matrices are the BCH 
formulas in SL(2,C), which can be obtained from the rephasing transformations in SU(2), by 
analytic continuation. We can represent this in a commutative diagram: 

BCH SU(2) 



N 2 



Analytic 
continuation 



BCH SL(2,C) 



Analytic 
continuation 



(37) 



N 2 



We now turn to the case of three flavors. A general symmetric SU(3) matrix (with five 
parameters) can be parametrized in either of two ways 
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vw diag v T , 
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The order of the non-commuting matrix products are chosen so that W\ corresponds to the 
usual mass matrix parametrization, while W 2 is most convenient for rephasing considerations. 
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The BCH formulas yield analytic relations between the two sets of parameters (r/j,e,) and 
A rephasing transformation would change 5j — > 5j, and the corresponding transforma- 
tion on (r]i, 6j) can be calculated as in Eqs.( p7H30| ). Also, the functions £j in terms of (r/j, ej) are 
rephasing invariants. The functions for (r]i,ej) obtained would have provided explicit solutions 
to the RGE, as in Eq.(^). Unfortunately, owing to the complexity of the SU(3) algebra, so far 
we are unable to solve for these functions explicitly. 

When we let all parameters assume complex values, the matrix W turns into a symmetric 
mass matrix. A rescaling (renormalization) transformation corresponds to 5j — > 6j +iAj. Thus, 
the rescaling (RGE) invariants which are also rephasing invariants are precisely (^1,^4,^6)) as 
complex functions of (rji , ej ) . Although we can not obtain these functions explicitly, we can 
obtain the RGE invariants using the matrix elements as the variables. As before, we can arrive 
at the results by first studying the rephasing invariants in SU(3). 

Consider rephasing transformations on a general SU(3) matrix (with eight parameters), 
which can be written in the form 

y = Pe~ ier ^ 7 e~ ie5 ^ 5 e~ ie3 ^ 3 e~ ie2 ^ 2 P' (43) 



where P is given in Eq.(42) and P' is obtained from P by the substitution Sj — > 6j. This is 
precisely the parametrization for the CKM matrix and (e 2 , £5, 63, ej) are rephasing invariants. 
In terms of the matrix elements of V, Vu, rephasing transformation gives Vjj — > e l ^ Sl+s J^Vjj. 
Thus, the rephasing invariants are |V/-j| 2 . When we let the parameters be complex, V analyti- 
cally continues into M, the mass matrix. The rephasing and rescaling invariants are now given 
by MfjMjj. 

Besides |V/j| 2 , another familiar form of the rephasing invariant ||| is given by Jijkl = 
VijVklVilVrj- When we impose the condition detV = 1, only relative rephasings, but no 
overall phases, are admitted. The rephasing invariant takes a simpler form 

h = e IJK ep j> k'Vii'VjjiVkk' . (44) 

There are six different ways to arrange the indices so that we may label / by s, which is an 
element of the permutation group S3, 

denoting (J I' , J -► J' ', K -► K'). 
Note that 

J2ls = detV = 1. (45) 

s 

Also, I s has a simple relation to the familiar rephasing invariant Jijkl- For a unitary V with 
det = 1, its minors are just the complex conjugated elements. For instance, 

V u = V 2 * 2 V 3 * 3 - V* 3 V£ 2 , or V n V 22 V 33 = |V2 2 | 2 [^ 33 | 2 - V 22 V 33 V 2 * 3 V^. (46) 




7 



When we let the indices take on different values, it is easy to show that all of the products I s 
are similarly related to Jijkl-, with Im/5 = —hn.{JijKL)i independent of the indices. 

The analytic continuation of the rephasing invariants I s turns them into rescaling and rephas- 
ing invariants for the mass matrices. 



Js = z-ijk &LMN Mil Mjm Mkn ■ 



(47) 



4 RGE Invariants 



As we discussed in the previous section, the physical RGE invariants are the rephasing and 
rescaling invariants of the mass matrix. There are two equivalent forms of these invariants. 
1) MjjMu, 2) 6jjkglmnMjlMjmMkn • When we convert these into physical variables, it 
turns out that the former is more convenient, which will be presented in the following. 
Let us define 

I T j = MjjMjj. (48) 
These invariants are not independent, since Mjj = Mjj, and 



/ j 

So there are altogether three independent (complex) invariants, which we can take to be 



(49) 



h =1 
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hh = I 
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I13 and h = I- 



23- 
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Explicitly, we find 
h 
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Here we have used the notations in Eqs.(lC-[ll)- In addition, 
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(51) 
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Ar 2 



^r 3 



X2 = e- 2 ^ 



I (e^ - e~ 2 ^) , (53) 

\{xi~e^), (54) 

( c 2 2 ) e ^+4) e " 2% )' (55) 

(s 2 2) e- 2 ^ + c 2 2) e-^) . (56) 



These invariants show that the physical parameters are intricately correlated during the RGE 
evolution. In general, we can not single out one or two variables which evolve independently of 
the others. However, in limited regions when certain conditions are satisfied, we do get simplified 
relations between a subset of the parameters. We will highlight some of these relations in the 
following. 

1) Real Mass matrix. 

In this case, all physical phases vanish, e 3 = 7^ = 0. The RGE invariants I\ 23 are all real 



so that renormalization does not generate any physical phases, as expected from Eq.(19). 

2) Two flavor solutions. 

The three flavor problem reduces to that of two flavor under certain conditions. This happens 
when two of the three mixing angles vanish. Thus, if = = 0, we find that I\ = J 2 = 
s^) sinh 2 ^ —rj 2 ),l3 = 0. Note that the condition S( 5 ) = S( 7 ) = is RGE stable. Similarly, if 
s (2) = s (7) = 0) or if s (2) = s (5) = 0, the result are genuine two flavor solutions. 

3) In regions when one angle is small. 

It may happen that, in a certain range of t, one of the mixing angles can be small. In general, 
such conditions can only be fulfilled in a limited region. There will then be approximate invariant 
combinations from a reduced set of parameters. For instance, in a region where S( 2 ) ~~ * 0> we 
find (Tfr = rj 1 - ie 3 ,% = r/ 2 + ie 3 ), 

J i -» s 2(5) sinh2 (^l -Va) 

h -> -C 2 ( 7 )S2(5) Smn2 (% -V3) ( 57 ) 



T 2 



, sinh 2 ( % - Vi) ~ C( 5 ) sinh2 ( r ?2 - V3) + s (5) c (5) s[nh2 (Vi ~ V3) 



This means that, in the limit S( 2 ) — * 0; ^he (1 — 3) sector behaves like a two flavor problem. 
At the same time, there is a correlation between the mixing angle ej and the CP phase, €3. 
However, if €3 = 0, then the approximation sr 2 ) — ► is consistent only if ej — > 0. Similar 
conclusions can be reached for the case s/q\ — > and sr?) — > 0. In addition, the cases when two 
masses are nearly degenerate, or when one mass value dominates, can also be analyzed along 
these lines. 
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5 Conclusion 



In this paper, we have studied the properties of the three flavor neutrino mass matrix under 
RGE evolution. Unlike the two flavor problem, where we obtained exact analytic solutions, the 
algebra of the 3x3 matrices is formidable, and we are only able to find three (complex) RGE 
invariants which correlate the evolutions of the many physical parameters. 

The RGE evolution of a (symmetric) mass matrix with detM = 1, for FCNC conserving 
theories (including the SM and MSSM), is given in the form M — > e^Me®, Q = real diagonal and 
trQ = 0. That is, renormalization amounts to a relative rescaling between the different flavors. 
At the same time, M is subject to arbitrary rephasing transformations, which correspond to 
taking Q to be an arbitrary pure imaginary matrix. The combined rescaling and rephasing 
transformation is thus given by M — » e^Me®, Q = complex. Since Q contains only two 
(complex) variables, we expect RGE invariants formed from the many physical variables in M . 
By means of an analytic continuation, these considerations are the same as those in obtaining 
rephasing invariants of the CKM matrix. We can thus write down three (complex) RGE and 
rephasing invariants explicitly. Our arguments also make it clear that these invariants are 
universal, independent of any specific model used to calculate the RGE evolution. 

Since exact solutions for the three flavor problem are not available, a number of approximate 
solutions have been considered in the literature B. The RGE invariants can be used to check 
the consistency of these approximations and to suggest viable new ones. The structure of 
these invariants also shows that, while the two flavor approximation is natural in a number of 
situations, their validity can only be established for a limited range of t. For large t, when 
the parameters also vary considerably, the two flavor approximation is viable only under very 
stringent conditions. 

With minor changes, most of the arguments in this work can be adapted to the study of 
quark mass matrices. In fact, pi it was shown that the infrared fixed point for two flavor 
RGE evolution corresponds to /3 — > Q,mi,/m\ — > oo. The approach to the fixed point, however, 
is governed by Eq. (|36|) , giving (i\frw^rn\ — > constant. This suggests that the quark mass 
matrices are the results of large RGE evolution, and that the well-known empirical relations 



between mixing angle and mass rations, O^j ^ *J Tn^j Trij , may have a dynamical origin. We plan 
to apply our analysis to a detailed study of the quark sector in the future. 
This work is supported in part by DOE grant No. DE-FG02-91ER40681. 
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